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We study the dynamics in phases with vector condensates of gluons (gluonic phases) in dense 
two-flavor quark matter. These phases yield an example of dynamics in which the Higgs mechanism 
is provided by condensates of gauge (or gauge plus scalar) fields. Because vacuum expectation 
values of spatial components of vector fields break the rotational symmetry, it is natural to have 
a spontaneous breakdown both of external and internal symmetries in this case. In particular, by 
£ — ' using the Ginzburg-Landau approach, we establish the existence of a gluonic phase with both the 

rotational symmetry and the electromagnetic U(l) being spontaneously broken. In other words, this 
phase describes an anisotropic medium in which the color and electric superconductivities coexist. 
£NJ ■ It is shown that this phase corresponds to a minimum of the Ginzburg-Landau potential and, unlike 

» I ' the two-flavor superconducting (2SC) phase, it does not suffer from the chromomagnetic instability. 

f"S |. The dual (confinement) description of its dynamics is developed and it is shown that there are 

• light exotic vector hadrons in the spectrum, some of which condense. Because most of the initial 

symmetries in this system are spontaneously broken, its dynamics is very rich. 
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I. INTRODUCTION 



It is expected that at sufficiently high baryon density, cold quark matter should be in a color superconducting 
, state p], 0, Q. On the other hand, it was suggested long ago that quark matter might exist inside the central region 
of compact stars This is one of the main reasons why the dynamics of the color superconductivity has been 
intensively studied (for reviews, see Ref. 1). 

Bulk matter in compact stars should be in /3-equilibrium, providing by weak interactions, and be electrically and 
color neutral. The electric and color neutrality conditions play a crucial role in the dynamics of quark pairing 
[1, 0, H 0, El, EL E2| • Also, in the dense quark matter, the strange quark mass cannot be neglected. These factors 
pH j lead to a mismatch Sfi between the Fermi momenta of the pairing quarks. 

As was revealed in Ref. [l3[, the gapped (2SC) and gapless (g2SC) two-flavor color superconducting phases suffer 
from a chromomagnetic instability connected with the presence of imaginary Meissner masses of gluons. While the 
8th gluon has an imaginary Meissner mass only in the g2SC phase, with the diquark gap A < Sfi (an intermediate 
coupling regime), the chromomagnetic instability for the 4- 7th gluons appears also in a strong coupling regime, with 
5(1 < A < \/25u. Later a chromomagnetic instability was also found in the three-flavor gapless color- flavor locked 
(gCFL) phase 0, El, Ell • 

Meissner and Debye masses are screening (and not pole) ones. It has been recently revealed in Ref. [l7| that the 
chromomagnetic instabilities in the 4- 7th and 8th gluonic channels correspond to two very different tachyonic spectra 
of plasmons. It is noticeable that while (unlike the Meissner mass) the (screening) Debye mass for an electric mode 
remains real for all values of dju, both in the 2SC and g2SC phases [l3|, the tachyonic plasmons occur both for the 
magnetic and electric modes [17] . The latter is important since it clearly shows that this instability is connected 
with vectorlike excitations: Recall that two magnetic modes correspond to two transverse components of a plasmon, 
and one electric mode corresponds to its longitudinal component. This form of the plasmon spectrum leads to the 
unequivocal conclusion about the existence of vector condensates of gluons in the ground state of two flavor quark 



'Electronic address: gorbar@bitp.kiev.ua 

^Electronic address: michioh@eken.phys.nagoya-u.ac.jp 

tElcctronic address: vmiransk@uwo.ca. On leave from Bogolyubov Institute for Theoretical Physics, 03143, Kiev, Ukraine 



2 



matter with A < y/25(i, thus supporting the scenario with gluon condensates (gluonic phase) proposed in letter [i~8| . 
While the analysis in Ref. [l8| was done only in the vicinity of the critical point 5a ~ A/y/2, a numerical analysis of 
the gluonic phase far away of the scaling region was considered in Refs. 120I w\ . 

At intermediate energy scales of the order of the diquark condensate A ~ 0(5OMeV), the analysis of QCD dynamics 
is very hard. Hence the phenomenological Nambu-Jona-Lasinio (NJL) model plays a prominent role in the analysis 
in dense quark matter (for recent extensive studies of dense QCD in this approach, see Refs. [UG3). The NJL 
model is usually regarded as a low-energy effective theory in which massive gluons are integrated out. The situation 
with dense quark matter is however quite different from that in the vacuum QCD. We will introduce gluonic degrees 
of freedom into the NJL model because in the 2SC/g2SC phase the gluons of the unbroken SU(2) C subgroup of the 
color SU(3) C are left as massless, and, under certain conditions considered below, some other gluons can be also very 
light. This yields the gauged NJL model. 

Also, because of the presence of matter, the running of the QCD coupling constant dramatically changes. As 
was shown in Ref. [13], the confinement scale Aq CD in the 2SC phase is essentially smaller than the typical scale of 
the diquark condensate, Aq CD < O(lOMeV) or even much smaller. This justifies introducing the gluonic degrees of 
freedom into the NJL model at the energy scale around and below A (and larger than Aq CD ), as was done in Ref. [l8l ]. 
At such scales, the dynamics in the gluonic phase corresponds to the Higgs picture. 

As we will discuss in this paper, the description of the infrared dynamics in the gluonic phase depends on the relation 
between two scales: Aq CD and the value of of the vector gluon condensates. If the latter is larger than Aq CD , then 
the Higgs picture is appropriate even in the infrared region, similarly as it happens in the electroweak theory. Indeed, 
when the SU(2) C gauge symmetry becomes completely broken by the dynamics with a characteristic scale being 
essentially larger than Aq CD , the strong coupling dynamics presented in the 2SC solution at the scale of order Aq CD 
is washed out. However, if this characteristic scale is < Aq CD , the confinement picture should be used. Although 
such a dynamics is not under control, the structure of the global symmetry should be the same both in the Higgs and 
confinement phases, if the Higgs fields (both vector and scalar ones) are assigned to the fundamental representation of 
the gauge group. In particular, the global charges of hadrons in the confinement picture should correspond to those of 
fundamental fields in the Higgs picture. This is one of the manifestations of the complementarity principle (23l. l24l. f25j] . 
Using this approach, it will be shown the existence of exotic hadrons in the gluonic phase. Moreover, we will see that 
dynamics with vector gluon condensates in the Higgs phase correspond to dynamics with condensates of exotic vector 
mesons in the confinement one. 

From the viewpoint of quantum field theory, it is quite natural to expect the existence of vector condensates of 
gluons in order to remove such instabilities as the chromomagnetic one or those connected with tachyonic plasmons. 
In Ref. homogeneous, i.e., independent of spatial coordinates, vector gluon condensates were considered. At 

the same time, because the condensates of spatial components of gluon fields break the rotational symmetry, such 
condensates are anisotropic. On the other hand, inhomogeneous condensates for diquark fields were studied in Refs. 
[li,[23,[Mj2lj33,|3l|. For example, the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase [H in QCD was considered 
m Refs. (illf - In this phase a diquark condensate is inhomogeneous. Some solutions with inhomogeneous diquark 



condensates, including the single plane- wave LOFF phase, can be considered as a special case of homogeneous vector 
condensates of gluons with the field strength being zero [H], [H, [34|. However, when the latter is not zero, as in the 
case of the gluonic phase [l8|, homogeneous vector gluon condensates cannot be traded for inhomogeneous diquark 
ones. Moreover, as was suggested in Ref. [l7|, inhomogeneous vector condensates of gluons could exist in the g2SC 
and gCFL phases. 

The paper is organized as follows. In Sec. [TT1 a renormalizable model with condensates of gauge fields is considered. 
This essentially soluble model yields a proof that such dynamics exist indeed. In Sec. IIII1 the gauged NJL model 
is described. Because the diquark gap A breaks the initial color SU(3) C symmetry to the SU{2) C one, it is useful 
to decompose the fields in the gauged NJL model with respect to the SU(2) C subgroup. This decomposition is 
considered in Sec. IIVI Section [V] is the central in this paper. In its four subsections, the dynamics in the gluonic 
phase is described in detail. One cannot exclude that besides this phase, other phases with vector condensates of 
gluons may exist in the gauged NJL model. In Sec. |VH we classify possible sets of homogeneous gluon condensates and 
typical symmetry breaking patterns in the corresponding phases in this model. We also describe the sets of operators 
relevant for constructing the Ginzburg-Landau effective theories for these dynamics. In particular, as a consequence 
of this analysis, it is shown that the ansatz for gluon condensates used for the gluonic phase [HI is self-consistent . In 
Sec. IVII[ the main results of the paper are summarized. In Appendices IA1 and IB1 , some useful formulas and relations 
are derived. 
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II. RENORMALIZABLE MODEL FOR DYNAMICS WITH VECTOR CONDENSATES 

Since a dynamics with vector condensates is a rather new "territory" , it would be important to have an essentially 
soluble model which would play the same role for such a dynamics as the linear <r models play for the conventional 
dynamics with spontaneous symmetry breaking with condensates of scalar fields. Fortunately, such a model exists: it 
is the gauged linear SU(2)l x U(1)y er-model (without fermions) with a chemical potential for hypercharge Y [3||. 1 

Let us briefly describe this model. It will be very useful for better understanding the dynamics in the gluonic phase. 
Its Lagrangian density reads (we use the metric g^ v = diag(l, —1, —1, —1)): 

C = -1fWfM«0 - If^F^W + \{D V - i^ Y 8 v0 )^{D v - i^yS" )^ - m 2 $^$ _ A($t$)2 ; (!) 

where the covariant derivative = <9 M — igA^ — (ig'/2)B fl , $ is a complex doublet field <f> T = (ip + ,ipo), and the 

chemical potential /iy is provided by external conditions (to be specific, we take fi Y > 0). Here = A^ ] T a /2 are 

SU{2) L gauge fields (r Q are three Pauli matrices) and the field strength F$ = d^A^ - d v A^ + ge abc A^ Ai c) . 

is a Uy(l) gauge field with the field strength Ffi? = d^B v — d u B^. The hypercharge of the doublet <E> equals +1. 
This model has the same structure as the electroweak theory without fermions and with the chemical potential for 
hypercharge Y. Note that the terms with the chemical potential are SU(2)l x U(l)y (and not SU(2)l x SU(2)r) 
symmetric. This follows from the fact that the hypercharge generator Y is Y = 21^ where /| is the third component 
of the right handed isospin generator. Henceforth we will omit the subscripts L and R, allowing various interpretations 
oftheSC/(2). 

The model is renormalizable and for small coupling constants g, g and A, the tree approximation is reliable there. 
Because the chemical potential explicitly breaks the Lorentz symmetry, the symmetry of the model is SU(2) x U (l)y x 
50(3) ro t- As was shown in Ref. [35[, for sufficiently large values of the chemical potential /zy, the condensates of 
both the scalar doublet $ and the gauge field A^ occur. The ground state solution is given by 

l(^)| 2 = ^-|, (4 3) ) = ^, <<f T ) = (0,,o), (2) 

where 



VCg 2 + 64AK - 8(8A - ff 2 )m 2 - 3, 9M y 
V ° = ^(8A-, 2 ) ' (3) 



= -^(A^ ± iA^), $ 3 = ((p + , (pa), and the vacuum expectation values of all other fields are equal to zero [39|. 

It is clear that this solution implies that the initial symmetry SU(2) x U(1)y X .SO^rot is spontaneously broken 
down to SO(2) rot . In particular, the electromagnetic U(l) em , with electric charge Q em — I3 + Y/2, is spontaneously 
broken by the condensate of W bosons, i.e., electric superconductivity takes place in this medium. 2 

Because the dynamics in this model is under control for small g, g' and A, the model provides a proof that the 
dynamics with vector condensate is a real thing. Moreover, this dynamics is quite rich. In particular, as was shown in 
Ref. [4j|, there are three types of topologically stable vortices in model ([T]), which are connected either with photon 
field or hypercharge gauge field, or with both of them. As we will see below, the dynamics in the gluonic phase 
strikingly resembles the dynamics in this toy model being however much more complicated. 



1 Ungauged linear 51/(2)^ X U(l)y o"-model with a chemical potential for hypercharge 36] is a toy model for the description of the 
dynamics of the kaon condensate in high density QCD |37l |. In particular, it realizes the phenomenon with abnormal number of Nambu- 
Goldstone (NG) bosons |36fl , when spontaneous breakdown of continuous symmetries leads to a lesser number of NG bosons than that 
required by the Goldstone theorem (for a recent discussion of this model, see Ref. [38j'). 

2 Note that because the U(l)y symmetry is local, for a nonzero chemical potential fiy one should introduce a source term BqJq in 
Lagrangian density JTJ in order to make the system neutral with respect to hypercharge Y . This is necessary since otherwise in such a 
system thermodynamic equilibrium could not be established. The value of the background hypercharge density Jo (representing very 
heavy particles) is determined from the requirement that Bo = is a solution of the equation of motion for Bo (the Gauss's law) [35, 40]. 
There exists an alternative description of this dynamics in which a background hypercharge density Jo is considered as a free parameter 
and fiy is taken to be zero. Then Gauss's law will define the vacuum expectation value (-Bo). It is not difficult to check that these two 

approaches are equivalent if the chemical potential fiy in the first approach is taken to be equal to the value ^r(Bo) from the second 
one. 
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III. GAUGED NJL MODEL 



We study dense two-flavor quark matter in /^-equilibrium. For our purpose, it is convenient to use a phenomenolog- 
ical NJL model with gluons, the gauged NJL model. As was already pointed out in the Introduction, this approach 
corresponds to the Higgs picture. The confinement picture will be considered in Subsec fVDl 

For simplicity, the current quark masses and the (V>V0 2 interaction in chiral channels will be neglected. Then the 
Lagrangian density is given by 



4 ' 



(4) 



where 



D„ = ft, - igA^T a , 



pO) = ft /i( a ) 



(5) 



Here A^f 1 are gluon fields, T a are the generators of St/ (3) in the fundamental representation, and f abc are the structure 
constants of St/(3). The spinor field i/j = ipi a has the flavor (i = u, d) and color (a = r, g, b) indices (ip ur is an up-red 
quark, etc.). In Eq. (e) = e lJ , (e ud = +1), and (e a ) = t al31 , (e rgb = +1) are the totally antisymmetric tensors in 
the flavor and color spaces, respectively. Note that we do not write explicitly the free electron term, and the photon 
field is not introduced in the model. 

In the /3-equilibrium, the chemical potential matrix /2q for up and down quarks is: 



(6) 



where 1 = l c <g> 1/, Q em = l c (g) diag(2/3, — 1/3) /, Q$ = diag(l/3, 1/3, — 2/3) c (S> 1/, and Me and (i$ are the quark, 
electron and color chemical potentials, respectively (the baryon chemical potential /is is fis = 3/i). Here the subscripts 
c and / mean that the corresponding matrices act in the color and flavor spaces, respectively. [Henceforth we will 
not show explicitly the unit matrices 1, l c , and 1/.] Note that the status of the chemical potential /ig is somewhat 
different from that of fi and fi e (in the absence of a photon field). The point is that the color neutrality condition is 



nothing else as the Gauss's law for the gluon field A y and /x 8 is expressed through A K q' as 



.(8) 



M8 



(7) 



As was shown in Ref. [42|, us is nonzero in the 2SC/g2SC phase. Although the color chemical potential is not an 
independent quantity in the gauged NJL model, we will keep the notation fig, in order to exhibit a special role of the 

field A, 



(8) 



Eq. ^ implies that the total chemical potentials for different quarks in the 2SC/g2SC phase are 



Mdr 



Mug 
Hdg 



: [1 
M 



Sfi, 



Hub 



M 
M 



Ms, 
Ms, 



(8a) 
(8b) 



with 



fx = fi- 



Let us now introduce the diquark field 

|$«|2 



4G A 



Me 
G 



Ms 
3 : 



Me 

2 



iip ee a ~f5ip. Then one can rewrite the Lagrangian density as 



4 



(9) 



(10) 



Without loss of generality, the diquark condensate in the 2SC/g2SC phase can be chosen along the anti-blue 
direction: ($ r ) = 0, = 0, ($ b ) ^ . In correspondence with that, in the gauged NJL model, it will be 

convenient to use the following (partly unitary) gauge in the 2SC/g2SC phase [43j |: 



0, $ 9 = 0, $ b 



A 



(11) 
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with the field A being real. Then the gap A in the 2SC/g2SC phase is equal to the vacuum expectation value of A, 
A=(A>._ 

We now introduce the Nambu-Gor'kov spinor, 



The inverse propagator of \& with the field A in the 2SC/g2SC phase is given by 

S~\P) = 
with 

[G+]-\P) = (po + p - Spr 3 - M8 1 6 ) 7 ° - 7 • p, 



[Gtr 1 a- 

A+ [Go]- 1 



and 



[G ]" (P) = (po - A + Sun + M&h -1-p, 



A- = -ise b j 5 A, A+ ee 7 °(A-) t 7° = -fee%sA, 



where P M = (po,p) is an energy-momentum four vector, t 3 ee diag(l, —1)/ and lb = diag(0,0, l) c . 
The propagator of \& is given by 

S(P) = (^G- 



with 



G ± = {[G t }- 1 -AT G TA ± }-\ H± 
The structure of G^ and was determined in the second paper in Ref. [l3| 



G f A ± G ± . 



G ± (P) ee diag(G±, G±,G±) C , 



where 



gup) 



Gt(P) 



(po T SfiT 3 ) - E ± 
(poT5pr 3 f-(E^f 



(poT6ht 3 ) + ET 



(po T 5pr 3 f - (Elf 



°A+ + 



(po T Spr 3 T Pa) + E ± p (p =f Spr 3 T Ms) - E^ 



-t°A- 



and 



with 



3*(P)EEe°( _I ± "la 



s± 



S±(P) = -<A 
S±(P) = -<A 



(po±^) 2 ~(^) 275 P 



A.+ + 



(poTW-(^) 275 p ' (p„ T <M 2 - 



A.+ + 



( P0 ±^) 2 -(^) 275Ap 
1 

2^ A P 



Here 



E ± = \p\±fi, E± = y/(E±Y + &, A±ee1(i± 7 



IpI J' 



(12) 

(13) 

(14) 
(15) 

(16) 

(17) 

(18) 
(19) 

(20) 
(21) 

(22) 

(23) 
(24) 

(25) 
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and while G^ and Gf are 8x8 matrices in the flavor-spinor space, the 4x4 matrices Sf 2 and 3^ act only in the 
spinor space. 

The generalization of expression (|13jl for the inverse propagator of '5 both with the scalar diquark field A and the 
vector fields in a gluonic phase (with nonzero gluonic condensates (A^) ) is straightforward: 

= ( [G |f ' [4,-. ) . M 

with 

[G+ g ]-\P) = (po + fioh°-i-P + g4 (a) T a , (27) 



[GoJ-^P) ee (po Ao)7° - 7 • P ~ gA^(T a f- (28) 
Integrating out fermion fields, we obtain the potential including both gluon and diquark fields: 

A 2 2 i n j4p 

T/ — I g f i°2°3 f«ia4«5 /l(a2) 4(03) 4(04) M /t(as) v I TV In f90A 

l 4G A + 4 / ; y 2j i{2*Y bs ■ { yj 

We will utilize the hard dense loop approximation, in which only the dominant one-loop quark contribution is taken 
into account, while the contribution of gluon loops is neglected. On the other hand, we keep the tree contribution 
of gluons in the effective potential (f2l))) . This is because we want to compare this contribution with that of hard 
dense loops in order to check the consistency of the hard dense loop approximation. The ground state of the system 
corresponds to the minimum of potential (|29[) and the question whether there exist gluon condensates (^4^) is a 
dynamical issue. 

The following remarks are in order. 

a) To study the ground state in the Higgs phase, it is convenient to use the unitary gauge. The important point 
is that in this gauge, all auxiliary (gauge dependent) degrees of freedom are removed. Therefore in this gauge the 

vacuum expectations values (VEVs) (Ap) of vector fields are well-defined physical quantities. The unitary gauge in 
the gluonic phase, which uses and extends the constraints presented in Eq. (JTTJ, will be described in Sec. [Vj 

b) As we will see below, in the gluonic phase, the time-component VEVs of the gluon fields other than the 8th 
one are also nonzero. Because of that, it will be convenient to rewrite effective potential (j2"9")l in a somewhat different 
form. Let us introduce the following matrix Ai g in the Nambu-Gor'kov space, 

MaT a 7° - gA a ■ 7 T a 

-Ma(T a ) T 7 ° + gA a -7 (T a ) 7 

where 



M g = S- x -S~ l = ^ 7 , -vr n ? ^ ST h (30) 



^ a = gAl, (S = 1,2,-.. ,7) (31) 
(/i 8 = ^-gA^ is included in 5" 1 ). Expanding now the logarithmic term in Eq. (f2"9")) , we find 

V = V A (A,p e ,m) + ^f aia2a3 f aiaia5 Al a ^Al a ^A^»A^ / — —Tr{SM g Y\ (32) 

4 n=l In J i{2n) 

where we defined the 2SC/g2SC part of the effective potential as 

V^,, & ) = §-- lJ^TrlnS-\ (33) 



The form (f32|) of the effective potential will be used in our analysis below. 



7 



IV. SU(2) C DECOMPOSITION 

Because the diquark condensate A = = (A) breaks the initial color SU(i) c symmetry down to the SU(2) C 
one, it is useful to decompose the initial fields with respect to the SU(2) C subgroup. In particular, the decomposition 
will help to calculate systematically the effective potential (15121) . 

The (anti-) fundamental and adjoint representations of SU(3) C are decomposed with respect to the SU(2) C as: 




Ipir 

1>ig 



®ip ib , (i = u,d), (34) 




3 = 291, i.e., = ( 35 ) 



and 

8 = 3®292®1, i.e., {Aft = (A™, A™ , A<®) 8 ^ 8 8 Af> , (a = 1, 2, ■ • ■ , 8). (36) 
Here we defined the complex doublets of the matter (with respect to the SU(2) C ) vector fields, 



j_ / 4 4) - *4 6) ) s = j_ f 4 4) + iA f 



(37) 



, / - - v ^ - iA)l> j y + 

Then we define the field strength for the SU (2) c gauge bosons, 

/« = - 9„A« + ge imn A^ A^ {i, m, n = 1, 2, 3) , (38) 

and the covariant derivative 

V^d^-igAf — (£=1,2,3). (39) 
It will be also useful to define the combinations 

A+ = + 4 2 '), -4; ^ - ^4 2 )). (40) 

Because of the presence of the electric chemical potential /i e in the model, the chiral SU{2) l^r symmetry is explicitly 
broken down to its J7(l)-part U(l) T z , Therefore the initial symmetry in the gauged NJL model is 

SU(3) e x U(l) em x [U{l) T s x U(1) T% ] X x SO(3) rot . (41) 

Since a photon filed was not included, the electromagnetic symmetry U(l) em is global. Note that the initial baryon 
charge is 

B= il/®l c = 2(Q em - 7 3 ), (42) 

where 1$ is a diagonal subgroup of the U{1) T 3 x U(1) T 3, i.e., ^3 = diag(l/2, —1/2)/. 
The diquark gap A breaks the initial symmetry (|41[) down to 

SU(2) C x U(l) cm x [U(l) T s x U(l) T s] x x SO(3) rot . (43) 
The unbroken U(l) era is connected with the new electric charge 
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where T 8 = ^Ajdiag(l, 1, — 2) c . The new unbroken baryon charge is 

£ = 2(Q em -J 3 ). (45) 

The charges for fermions, diquark fields and gluons are summarized in Tables HI and HT1 
The transformations of the fields under the gauge SU(2) C have the following form 



where 



<i> j I $'s j ' I $s ) • ( 47 ) 

0„ ^ = U<f> M , (48) 

4* -> ^ = C/^C/- 1 + -Ud^U- 1 , (49) 

-» = Uf^U- 1 , (50) 

p M p; t = uv^u~\ (5i) 



^=4^ y , f,u = f^\, 0?= 1,2,3) (52) 

and 

[/ = exp (i6 e [x)^j . (53) 

Under the new baryon symmetry connected with generator B (|45|) . blue quarks, the diquark fields <£ >,%9 and the 
field doublet transform as 

At -> V4 = e'Va, (i = u, d), * r ' 9 - $' r ' 9 = e* 9 <F' s , ^ -h. ^ = e~ i6 ^, (54) 

while other fields have zero baryon charge B (see Tables HI and IPTl ) 



V. GLUONIC PHASE 

Both the chromomagnetic [l3| and plasmon [l7| instabilities for the 4-7th gluons in the 2SC phase at A < y/2Sfi 
suggest a condensation of these gluons, i.e., the creation of a condensate of vector field (|37|) . Because the chromo- 
magnetic instability develops in the magnetic channel, it is naturally to expect that a spatial component of (j)^ has a 
VEV. In studying this condensate, one can use the freedom connected with symmetry (1431) in the 2SC/g2SC phase. 
By using the rotational symmetry SO(3) ro t, one can take (</> x ) ^ while ((p x ) = (<fi y ) = 0. And because of the SU(2) C 
symmetry, without loss of generality, we can choose (A^) ^ 0. This VEV breaks the SU(2) C down to nothing and 
the 50(3) rot down to the SO(2) mt . 

The following remarks arc in order, a) The complex doublet <j> z plays here the role of a Higgs field responsible 
for spontaneous breakdown of the SU(2) C . The situation is similar to that taking place in the electroweak theory. 
The essential difference however is that now the Higgs field is a spatial component of the vector field leading also to 
spontaneous breakdown of the rotational symmetry, b) In this paper, we will use the gauge in which 



^4=(o,<4 6) > + <4 6) ), 



(55) 
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TABLE I: The quantum numbers of the up and down quarks and the diquark fields. In the unitary gauge the diquark fields 
and Im& b are absorbed into the longitudinal modes of the corresponding giuons. 
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TABLE II: The quantum numbers of giuons. 



where the real field a z describes quantum fluctuations. This constraint together with that in Eq. (fTTj) constitute 
the unitary gauge: all auxiliary (gauge dependent) degrees of freedom are now removed. 

A gluonic phase with such a condensate was described in letter [l8|]. Since the most of the initial symmetries 
are broken in this phase, its dynamics is rich and complicated. On the other hand, because of space shortage, the 
description of this phase in letter [l8| was rather brief. In this paper, we will present both a detailed description of its 
dynamics and present a general analysis of a possibility of the existence of other phases with vector gluon condensates 
in dense QCD. 
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A. Symmetry breaking structure and Ginzburg-Landau effective potential 



Let us describe symmetry breaking structure in the gluonic phase. With a broken SU(2) C , the SU(2) C gluons could 
have VEVs. A similar situation takes place in the gauged a- model with a chemical potential for hypercharge described 
in Sec. [IT] above: the gauge symmetry SU(2)l is broken there. Motivating by that model, we assume 

(4 1 )), (4 3) )^0. (56) 

As will be shown below, a solution with these vector condensates exists in the model indeed. 

The symmetry in the 2SC/g2SC phase is that presented in Eq. {43]). The VEV (A^) breaks the SU{2) C but a 
linear combination of the generator T 3 from the SU(2) C and Q e m, 

Qcm — Ocm T Qcm ^ ' (^^) 

determines the unbroken C/(l) em - The new baryon charge is B — 2(Q cm — I 3 ) [the charges Q om and B for quarks, 

diquarks, and gluons are shown in Tables HI and ITT] . However, because T 1 does not commute with T 3 , the VEV {Ajp) 

breaks U em (l). The U(l) symmetry connected with the baryon charge B is also broken. 
After all, we have: 

[SU(3) c ] local x [U(l) cm x U(l) Tl x C/(l) T 3] globa i x SO(3) rot 

[SU(2) c ] local x [U(l) em x U(l) Tl x C/(l) T 3] global x 50(3) rot (58) 

^ [U(l) cm x U(l) T s x C/(l) T 3] g i oba i x SO(2) Tot (59) 

{A ^ ] [U(l) T s x C/(l) r 3] global x SO(2) Iot . (60) 

Thus, this system describes an anisotropic medium in which both the color and electric superconductivities coexist. 

Let us apply the Ginzburg-Landau (GL) approach to this system near the critical point 5fi ~ Aj\[2. The two 
point function of gluons can be calculated from Lagrangian density (fT0|) . While in Ref. [HI the Debye and Meissner 
screening masses of the gluons in the 2SC phase were calculated, the pole masses of the corresponding light plasmons 

(with masses \M\ <C were analyzed in Ref. [T^. For the gluons of the unbroken SU(2) C , i.e., A^\ A^\ and A^\ 
both Debye and Meissner masses vanish in the region 5fx < A and there are no light plasmons in these channels. For 

imat 



the gluons A^ 1 the Meissner mass is approximately 



6tt 2 



6(i < A. (61) 



Thus, near the critical point 5/j, = A/y/2, the Meissner mass for A^ ^ is very small. As <5/i exceeds the value 
A/y/2, m 2 M 4 becomes negative, thus signalizing a chromomagnetic instability in the 2SC solution. The pole masses 
of the light plasmons for the magnetic and electric modes have similar behavior in these channels. On the other hand, 
around the critical point d/u, = A/y/2, the SU(2) C singlet gluon A^ is heavy. Actually, it is heavy in the whole region 
Sfi < A. This fact allows us to pick up the gluons A^~^ as relevant light degrees of freedom in the low energy 
effective theory around the critical point Sfi = A/y/2. 

Because the SU(2) C is a gauge symmetry, the building blocks of the GL effective action are 

0o, fa, T> , Vj, f j, fjk, (62) 

where the indices j and k represent spatial components. The SU(2) C and 50(3) ro t symmetries dictate the form of 
the general GL effective potential, which is made from these building blocks and includes operators up to the mass 
dimension four, i.e., relevant and marginal ones. It will be convenient to introduce the following notations: 

B^gAf\ C^gAi}\ D = gAf . (63) 
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Then, as shown in detail in Subsec. IVI Bl below, the GL potential has the form 

V eS = V A + -M%B 2 + T DB DB 2 + X -X BC B 2 C 2 + X -X BD B 2 D 2 + ^X CD C 2 D 2 + \x B B\ (64) 

where V\ is the 2SC part of the effective potential (see Eq. (|33|) V Here, while the coefficients X B , Xbc, X B d, and 
Xcd are dimensionless, the dimension (in mass units) of the coefficient T BB in the triple vertex is one. Expanding 
the potential V (f3"2")l with respect to B, C, and D, we can determine these coefficients. 

Before realizing explicit calculations, we clarify the behavior of the effective potential ([64]) near the critical point. 
The stationary point of the effective potential (|64|) is given by the equations 

— J21 = B [ M% + X B B 2 + 2T DB D + X BC C 2 + X BD D 2 ] = 0, (65) 

dVcs 

dC 
dV cS 

3D 

and 



= C [ X BC B 2 + X CD D 2 ] = 0, (66) 
= T DB B 2 + X BD DB 2 + X CD C 2 D = 0, (67) 



dV eS dV cS dV cS 

^7 = °' ^7 = °' ^a =0 - (68) 

It will be convenient to present fj, e , fis, and A as 

Me = Me + £e, (69) 
M8 = (70) 

A = A + U, (71) 

where the bar-quantities are from the 2SC solution, with B = C = D = 0. Let us assume that the origin (bifurcation 
point) of the solution with nonzero B, C, and D corresponds to a second order phase transition (as will become clear 
below, this assumption is self-consistent). Under this assumption, the analysis of Eqs. (f6"S"]) - (f6"7| and (|6"5|) was done 
in Appendix [XJ Taking an infinitesimally small B near the critical point, it is shown there that 

~0(B 2 ) (72) 

and that when the 2SC solution becomes unstable (Mg < 0), a new solution occurs, if the parameters X B c and Xcd 
satisfy 

X BC > 0, X CD < (73) 
(in the next section, it will be shown that this constraint is satisfied indeed). The new solution is: 



" ''^5^1^ «-^">-ViS *—<^>=^ ™ 



where we neglected higher order terms in M^. It is important that, as shown in Appendix lAl the coefficients M|, 
Tdb, Abc, and Xcd in this nearcritical solution are expressed through the 2SC values A, j2 e and p,$. Note that in 
Eq. (174)) the convention B > and C > is chosen. 

Near the critical point M% = 0, the solution behaves as 



B sol oc — Mg, C sol cx ^J-M 2 , D sol oc —M B . (75) 

These scaling relations are quite remarkable. While the scaling relation for C is of engineering type, those for B and 
D are not (the origin of this is of course in the presence of the dimensional coefficient T BB in Eq. ([74]) ). Such a 
scaling behavior implies that the B A and B 2 D 2 terms in the effective potential are irrelevant near the critical point 
= 0. Omitting them, we arrive at the reduced effective potential: 

V cS = V A + \m 2 b B 2 + T DB DB 2 + l -X BC B 2 C 2 + ^X CD C 2 D 2 . (76) 
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Let us now turn to the 2SC part Va in Eq. (|76|) . As shown in Appendix IA1 the difference of Va in the new solution 
and that in the 2SC one is 

V A (A S °\ tf 1 , ^f) - F A (A, fl e ,fl s ) ~ 0(B 4 ). (77) 

On the other hand, as follows from Eqs. ([75]) and ([75]) , the difference Veg — Va ~ 0(-B 3 ). This fact and Eq. 
(|77p imply that in the leading approximation one can use the 2SC bar-quantities in calculating Va in the reduced 
potential. In other words, the effective potential can be decomposed into the "constant" 2SC part Va, with frozen 
fermion parameters, and the dynamical gluonic part: 

V- eff - Ve ff (A, p e ,/2 8 ; B, C, D) = V A {A, /2 e ,/2 8 ) + \m 2 b B 2 + T DB DB 2 + l\ BC B 2 C 2 + ^X CD C 2 D 2 . (78) 
Eq. (j74|) is the exact solution for the potential ([78]) and the energy density at the stationary point is 

V cff (A, fl e , j2 s ; B sol , C soh D sol ) = Va + ^/ 2 B 2 ol = V A - ) < ^ (79) 

Therefore the gluonic vacuum is more stable than the 2SC one. 

In the description of the dynamics with vector condensates, there is a subtlety connected with the derivation of a 
physical effective potential, whose minima correspond to stable or mctastable vacua. The point is that although the 
gauge symmetry is gone in the unitary gauge, the present theory still has constraints. In fact, it is a system with 
second-class constraints, similar to the theory of a free massive vector field described by the Proca Lagrangian (for 
a thorough discussion of systems with second-class constraints, see Sec. 2.3 in book [44[). In such theories, while the 
Lagrangian formalism can be used without introducing a gauge, the physical Hamiltonian is obtained by explicitly 
resolving the constraints. 

In our case, this implies that to obtain the physical effective potential, one has to impose the Gauss's law constraint 
on the conventional effective potential V c g (|T8[) . This constraint amounts to integrating out the time- like components 
A^\ In the present approximation, the latter can be done by using their equations of motion, which are reduced to 
Eq. (|57|) for D = gA^ in our case. Omitting the suppressed DB 2 -teim in this equation, we get 

T DB B 2 + \ CD C 2 D = 0. (80) 

It leads to the physical effective potential without the non-dynamical degree of freedom D: 

V c ^ s = V A + \M 2 B 2 + \\ BC B 2 C 2 - (81) 

It is easy to show that solution ([71)1 is a minimum by analyzing the curvature of V3c ya - Note that because of the 
constraint in Eq. (|73p . this potential is bounded from below. 

In the next section, we will calculate M B , T B b, A^Cj an d ^cd- In particular, it will be shown that constraint (|73[) 
is satisfied near the critical point. 



B. Dynamics in one- loop approximation 

In this subsection, we determine the GL effective potential (|78p in one-loop approximation and derive the dispersion 
relations for quarks in the gluonic phase. The 2SC Va part of the potential is known [101 ] . 



A 2 
4G A 

^2 



/' 4 



12tt 2 

-2 



(1 



l -A 2 
4 



12tt 2 12tt 2 3vr 2 
4A 2 A 2 



In— TT 7T 



A z 



2fl 2 



l -A 2 



(5 fx < A). 



(82) 



Here A is the ultraviolet cutoff in the NJL model and fi u b, Hdb, and jl are given in Eqs. © and ©. For clarity of the 
presentation, the bars in A, /i e and /xg were omitted [0(jl 2 /A 2 ) and C(A 2 /A 2 ) and higher terms are neglected in this 
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expression]. 3 Note that the color and electrical charge neutrality conditions in the 2SC solution yield [1C 



r 3 1 

10^5^' 



and 



{ff + 5fi 2 )fi s = - M A 2 In 



2A 



1 ) +/i(d> 2 + /i 2 )--/i 3 , 



(83) 



(84) 



which is consistent with the result of Ref. [13], Ms ~ 0(A 2 //i), in the case of (5/x = 0. The size of the diquark gap A 
is essentially determined by tuning the NJL coupling constant Ga and cutoff A. 

In Appendix [B] after straightforward but tedious calculations of relevant one- loop diagrams from the fermion trace 
in Eq. ([29]) . we find the following relations in the region 8fi < A: 



r,2 



^CD 



80tt 2 A 2 



-1 + 8 



5jl_ 
A 2 



Tdb — 



Ms 



18tt 2 A 2 : 



Ms M 2 
24tt 2 A 2 



-1 



A 4 



48tt 2 



-1 + 4 



Sfi 2 5^ 
A 2_ + 8 A4- 



(85) 
(86) 
(87) 



Here the tree contribution of gluons 



Va 



-Co 



L 77 " ; f 

2 °i °J 



(a) p (a) 



v^ 2 + v^ s2 



8ff 2 



r B 2 ,D 2 



2 r.2 



2<? ; 



:C 2 L> 



was also taken into account. As to the coefficient M§, its expression follows directly from Eqs. (l6"T]) and 



M% 



( -/i 2 + m 2 M 4 ) 



Ms ! M 
q 2 67T 2 



2V 
"A 2 " 



(89) 



We sec that the coefficient Acd in ([56]) is definitely negative. The parameter M|, which is expressed through the 
Meissner mass (|6"Tj) . is negative when 



6/i > <5/i c 
Relation ^ and Eq. © yield 



^Mcr = 



7T 



1 - 



3?r m| 

2a s /i 2 ' 



47T 



(90) 



M= 

and, at the critical point, we find from Eqs. (j83[) . (|84| and ([90"]) , (j9Tj) that /is is approximately 



(sol) 

Ms 



o ln 200A J 



12 



l n 222Ai_ 2 



M • 



For realistic values A = (1.5 — 2.0)/i and a s = 0.75 — 1.0, we obtain numerically 



3tt m|i 

2a s /2 2 'm8=M8 S ' 



0.03-0.1 



(91) 



(92) 



(93) 



3 For realistic values A = (1.5 — 2.0)/i that we use, while the contribution of the 0(A 2 /A 2 ) terms are parametrically suppressed, one can 
show that the contribution of the 0(fi 2 / 'A 2 ) terms is numerically suppressed. 



14 



This implies that the tree gluon contribution in Eq. ((90)) decreases the value of <5/x cr by 1.5%-5% in comparison to 
its value in the (non-gauged) NJL model. The smallness of this correction is in accordance with the dominance of 
hard-dense-loop diagrams. 

Let us now turn to the coefficient Xbc (|85|) . At the critical point 5(i — <5/i cr , it is [henceforth we will not show 
explicitly the superscript (sol) in fig,]: 



Xbc — 



1 



F j2 



80tt 2 A 2 



2« 2 ^ 



(94) 



Because the /^|//i 4 -term is negligibly small, we conclude that the coefficient Xbc 1S positive near the critical point. 
Thus, constraint (|73| is satisfied indeed. 

Utilizing Eqs. ([55]) - ([57)) and (|89[) in Eq. (f?4")) . one can obtain the solutions for B, C, and D in the near-critical 
region. Indeed, neglecting higher order terms in fi%/ fi in Eqs. (|85|) - ([87|) and ([89]) . we get the approximate relations 



M 



f' 2 



B 



1 - 



Sfi 2 



A 



9 



which lead us to the near-critical solution: 



A 



1 



CD 



Aira s 



1 

irT 2 ' 



Tj 



DB 



fl 



M8 



16tt 2 16tt 2 



o 27r 

3 + — 



(95) 



Sfi 2 - Sn 2 c 



16 yU 



Si4 



27 



^ f 3+ 2il 



(96) 



Csol 



(97) 



/ 3+ 2^ 

f ( as 



It is noticeable that this solution describes nonzero field strengths F^) which correspond to the presence of non- 
abelian constant chromoelectric-like condensates in the ground state: 



(4 2) ) 
(EP) 



~ C so \D so \ , 
9 

— B so i ( 2{i s — D so \ ) 
25 



(99) 
(100) 



We emphasize that while an abelian constant electric field in different media always leads to an instability, 4 
non-abelian constant chromoelectric fields do not in many cases: For a thorough discussion of the stability problem 
for constant SU(2) non-abelian fields in theories with zero baryon density, see Ref. (45|. On a technical side, this 
difference is connected with that while a vector potential corresponding to a constant abelian electric field depends 
on spatial and/or time coordinates, a constant non-abelian chromoelectric field is expressed through constant vector 
potentials, as takes place in our case, and therefore momentum and energy are good quantum numbers in the latter. 

In order to illustrate the stability issue in the gluonic phase, let us consider the dispersion relations for quarks there. 
Because the vacuum expectation values (|96[) - (198l) are small near the critical point and because red and green quarks 
are gapped in the 2SC phase, the dispersion relations for gapless blue up and down quarks are of the most interest. 



4 In metallic and superconducting media, such an instability is classical in its origin. In semiconductors and insulators, this instability is 
manifested in an creation of electron-hole pairs through a quantum tunneling process. 
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From Eq. ((26)) we find that up to the first order in B 2 they are 



rfl -ldl-„ | B sol 1 B 2 ol (p*y 1 2(|p1-/l)+M8 nnn 

P *- m -' h * + 4 2b1+/i8 + 5 AL_- 4 p 2 ^ 2 |pl +M8 + 5 ^ + A 2 - / ,I-2Mb1-A)J' (101) 



Pdb " b1 — + ^ 8 + ~ ~ ~fT ~ — ~ _0_ + a 2 - - 2M\p\ - A) ) ' (1 ° 2) 

The i? 2 -terms in Eqs. (|10ip and Q102p lead to non-spherical Fermi surfaces determined by the following equations: 



4 2 ^ + MB + ^^ 4 A 2 + /isA^e + m| 
N = Mdfc-^i X3 + - §2i l A2 blue down 104 

where we neglected higher order terms in B 2 and defined the angle 0, 

p 3 = \p\cos9. (105) 

The dispersion relations (|101[) and p02[) clearly show that there is no instability in the quark sector in this problem. 

As to bosonic degrees of freedom (gluons and composite bosons) , because it is very involved to derive their derivative 
terms from the fermion loop in the gluonic phase, this issue is beyond the scope of this paper. It is however noticeable 
that there are no instabilities for bosons in a phase with vector condensates in the gauged cr-model with a chemical 
potential for hypercharge [35j. Although that model is much simpler than the present one, its phase with vector 
condensates has many common features with the gluonic phase and this fact is encouraging. Note that among light 
collective excitations in the gluonic phase, there should be Nambu-Goldstone bosons connected with the spontaneous 

breakdown of the 50(3) ro t and the U(l) cm (in the presence of photon field, the latter will be absorbed in the electric 
(longitudinal) part of the field). 

We emphasize that these constant color condensates in the gluonic phase do not produce long range color forces 
acting on quasiparticles. This can be seen from the dispersion relations (jlOip and (|102[) for quarks in this model. 
They show that momentum and energy are conserved numbers. It would be of course impossible in the presence of 
long range forces. The role of these condensates is actually more dramatic: They change the structure of the ground 
state, making it anisotropic and (electrically) superconducting. Only in this sense, one can speak about a long range 
character of the condensates. 



C. Searching for other solutions 

Are there solutions of Eqs. ([65]) - (f67|) other than that found in the previous section and the trivial one with 
£>soi = C so i = D so \ = 0? We will address this question in this subsection. 

First of all, it is easy to see that there are two formally nontrivial solutions with B so \ = 0: a) C so i = 0, _D so i-arbitrary 
and b) D so i = 0, C so i-arbitrary. However, as follows from Eq. (JHSJ), both these solutions lead to zero SU(2) C field 

strength fyl ([38]) and therefore they are gauge equivalent to the 2SC solution without vector condensates. 

It is also easy to check that for D so \ — 0, Eqs. (f6"5|) - (f6"T)) lead to B so \ — 0, C so \ = 0. Therefore we conclude 
that both B and D have to be nonzero in a nontrivial solution. Physics underlying this conclusion is clear. When 
B so \ = g(Ag ) ^ 0, the color charge densities of red and green quarks are generally different, so that D so \ = g(A^) ^ 
is also required for color neutrality (recall that gA$ can be considered as a chemical potential ^3 related to the third 
component of the color isospin). 

There does exist one additional nontrivial solution with nonzero (E^ ) = (-F03 ) ( see Eq. (|100|l ): 

C sol = 0, D Bol = -^, B 2 ol = -^ + ^^. (106) 

*BD *B AB*BD 
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This solution corresponds to a phase in which while the gauge SU(2) C symmetry and the rotational SO(3) ro t are 

broken, the electromagnetic U(l) em is exact. While physics in this phase is quite interesting, there is the following 
problem in justifying its existence. As follows from Eq. (|95| . the coefficient Tdb does not approach zero at the 
critical point M B = 0. This and Eq. (|106[) imply that the values of B so \ and D so \ in this solution are also nonzero at 
the critical point. Therefore the solution corresponds to the first order phase transition and the GL approach is not 
appropriate in this case. Since the derivation of Eqs. (I65l) - (j67]) . (|78|) . and (f8Tj) was based on this approach, all these 
equations themselves will be modified in this phase. Therefore the question concerning the existence of solution (|106|) 
is open. It seems that a numerical analysis would be the only reliable way to answer it. 

In conclusion, the following two remarks are in order, a) Since both solutions (|96|) - (|98|) and (|106|) are cylindrically 
symmetric, it will be appropriate to call the phase corresponding to solution (|96p ~ (|98|) as a gluonic cylindrical phase 
I and the phase corresponding to solution (|106|) as a gluonic cylindrical phase II. b) The fact that any nontrivial 
solution should have both B and D to be nonzero follows from the presence of the triple vertex TdbDB 2 in the GL 
potential: for B ^ 0, this vertex inevitably leads to a nonzero D. One can call this a tadpole mechanism: For a given 
B, the diagram corresponding to such a vertex is a tadpole with a fermion loop (with B insertions), producing the 
coefficient Tdb, and with a tail being the field D. 

D. Confinement picture and exotic hadrons in gluonic phase 

In this subsection, we will describe some additional features of the gluonic phase. In particular, we will describe 
the confinement picture, which can be appropriate for the description of its dynamics in the infrared region, with the 
energy scale of order Aq CD < O(10 MeV) (or even much smaller) [22l |. As we will see, the dynamics in this dense 
medium includes light exotic vector mesons some of which can condense. 

The gluon condensates are mostly generated at energy scales between the confinement scale in the 2SC state A'q CD 
and the quark chemical potential, which is about 300-500 MeV. It is the same region where the chromomagnetic 
instability in the 2SC phase is created and where the hard dense loop approximation is (at least qualitatively) 
reliable. At such scales, gluons are still appropriate dynamical degrees of freedom and utilizing the Higgs approach 
with color condensates in a particular gauge is appropriate and consistent: It is a region of hard physics. Because the 
gluonic cylindrical phase I occurs as a result of a conventional second order phase transition, the gluon condensates 
are very small only in the immediate surroundings of the critical point Sfi = A/ a/2- Outside that region, their values 
should be of the order of the typical scale 5/j, ~ A ~ 50 — 100 MeV. As to the gluonic cylindrical phase II, because it 
is connected with a first order phase transition, the situation depends on whether it is a strong or a weak one. While 
for the former, one could expect that the gluon condensates are of order Sfi ~ A even at the nearest surrounding 
of the critical point [see Eq. (|106p ], for the latter, they could be of order ^'qc D there, modulo the question of the 
existence of this phase (see the discussion in the previous subsection) . 

The description of the dynamics in the gluonic phases in the infrared region depends on the value of the gluon 
condensates. If they are essentially larger than A'q CD , then the Higgs description is appropriate even in the infrared 
region, similarly as it happens in the electroweak theory. Indeed, when the SU(2) C gauge symmetry becomes com- 
pletely broken by the dynamics with a characteristic scale being essentially larger than A'q CD , the strong coupling 
dynamics presented in the 2SC solution at the scale of order Aq CD is washed out. In this regard, the gluonic phases 
with large vector condensates are similar to the color-flavor locked (CFL) phase with a large quark chemical potential 
fi, where the color condensates (although not vector ones) completely break the SU(3) C color gauge symmetry [46j. 

But what happens if the gluon condensates are small, < Aq CD ? This regime corresponds to the nearcritical 
dynamics in the gluonic cylindrical phase I and the confinement picture should be appropriate for the description of 
the infrared dynamics in this case. 

In order to answer this question, note the following. As one can see in Tables U and HH the electric charge Q cm and 

the baryon number B are integer both for gluons and quarks. Do they describe hadronic-like excitations? We believe 
that the answer is "yes" . The point is that in models with Higgs fields in the (anti-) fundamental representation of the 
gauge group, there is no phase transition between Higgs and confinement phases [231 124 , |25| and this is the case in the 
present model. Indeed, in the 2SC phase, the breakdown SU(3) C — > SU{2) C is triggered by the diquark condensate, 
which is assigned to the anti-fundamental representation of the SU(3) C , and gauge SU{2) C symmetry breaking occurs 
when the SU(2) r doublet vector field M develops the VEV. Because of that, we can apply the complementarity 
principle (23l . I24. |25| for the description of the dynamics in the gluonic phase with small condensates in the infrared 
region. What matters is the existence of the unitary gauge given by constraints (fTTj) and (|55|) . In this gauge all gauge 
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dependent degrees of freedom are removed. 

Due to the complementarity principle, the Higgs and confinement phases provide dual, and physically equivalent, 
descriptions of dynamics. In particular, they provide two complementary descriptions of a spontaneous breakdown 
of global symmetries, such as the rotational SO(3) and the electromagnetic £ 7(1) in the gluonic cylindrical phase I 
and the rotational SO(3) in the gluonic cylindrical phase II. Following Ref. [25j, we will consider the dual, gauge 
invariant, approach in this model and show that all the gluonic and quark fields can indeed be replaced by colorless 
composite ones. 

The flavor quantum numbers of these composite fields are described by the conventional electric and baryon charges 

Q em and B. They are integer and coincide with those the operators Q em and B yield for gluonic and quark fields. 
The composite fields in confinement picture should coincide with the corresponding fields in the Higgs picture in the 
unitary gauge in the classical approximation. 

The nonzero VEVs in the Higgs picture (common for these two gluonic phases) are: 



(*) = (0,0,A) T , fi 3 = D sol =g(A^), M8 



(4 



(8)\ 



B, 



sol 



(107) 



and thereby 6 



(iD z &*) = (0, B sol A/2, Of, g{F 7 0z ) = B sol Us - ^ ) , g(F 0z <*>*) = 



-i%^(2 M8 -/x 3 )(0,l,0f, (108) 



where 



p = p(a)rpa rpa 

r jJ,V — r fj, v 1 ± — 2 



(109) 



with A a 's being the SU(3) C Gell-Mann matrices. In the classical approximation, we replace the above fields by their 
VEVs. Then the following composite fields can be written in terms of the elementary fields: 



/4 4) - ^4 5) \ 



-—A 



A 



(6) 



iA 



(7) 



_ 4(8) 

\/3 J 



(110) 



( 4 1} 



^^S so iA 



iA 



(2) 



-A 



(3) 



V3A 



(8) 



2iA 



(7) 



(111) 



and 



/(2/is 



4 



iA 



(5) 



(2 M8 - M3 )(4 6) 



iA 



(7)> 



V 







/ 



(112) 



By using the above relations, we can construct composite fermions and bosons in confining picture. For example, 
the (up and down) blue quark fields can be rewritten as $> T i(>i with "0; = (ipi r ,ipig,'4'ib) T , (* = u,d). Note that 



5 Recall that we do not show explicitly the superscript (sol) in the chemical potentials fi^ and [ig. 

6 Recall that the diquark field <I> is an anti-triplet under the SU(3) C symmetry. 
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~ Re$ 6 











(F 0z **)t(F 0z **) ~Re$" 












TABLE III: Composite fermions and scalar fields in the confinement picture. 

in the classical approximation the composite fields & T ijJi yield Aipn,. The green quarks in the confinement picture 
can be described by (D z $*)^ipi — > ^B so iAipi g . By using the epsilon tensor, we can rewrite the red quarks as 
e a °"t'ip ia (D z $*)i3($*) 1 — > 7}B so \A 2 ipi r . For the diquark field, only the real part of the anti-blue one is a physical 
degree of freedom. The composite field is €>^$ — > 2ARe$ b . For other fields, see Tabic ITTT1 

Similarly, we can construct the vector composite fields. For example, we find (D Z &*Y (D^*) — ► \B so \A 2 (A^' — 
iA { J ] ) ~ 09, e «^(^**) Q ( J D z **) (**) 7 -> iB sol A 3 (4 4) - iA^ ] ) ~ ^, e^ 7 (F 2M **) Q (^**) /3 (**) 7 -> 

|Bg ol A 3 (^4^ - M^ 2) ) - A~, etc.. We summarize them in TableE] 

Some of vector mesons in Table ITVl are exotic because they carry baryon charge. (In vacuum QCD, mesons carry of 
course no baryon charge.) For example, the electric and baryon charges Q om and B of vector mesons corresponding 

to = Ap ± iA^ gluons are equal to ±1 and ±2, respectively. The origin of these exotic quantum numbers is 
connected with (anti-) diquarks, which are constituents of these mesons (see Table HVl) . Indeed, (anti-) diquarks are 
bosons carrying the baryon charge ±2/3 and therefore are exotic themselves. 

This feature has a dramatic consequence for the gluonic cylindrical phase I. Since in the Higgs description of this 

phase A( ±S) gluons are condensed (leading to the spontaneous U (l) e m breakdown), we conclude that in the confinement 
picture this corresponds to a condensation of exotic vector mesons. In this regard, it is appropriate to mention that 
some authors speculated about a possibility of a condensation of vector p mesons in dense baryon matter |47j. The 
dynamics in the gluonic phase yields a scenario even with a more unexpected condensation. 

VI. DYNAMICS WITH GLUON CONDENSATES: GENERAL ANALYSIS AND CLASSIFICATION 

In the previous sections, the dynamics with gluon condensates connected with the instability for the 4- 7th gluons 
in the 2SC phase was considered. The question is whether there are gluonic phases other than those two described 
in Sec. El In particular, an interesting issue is the dynamics of gluon condensates connected with the instability for 
the 8th gluon in the g2SC phase. Moreover, our consideration of the gluonic phase in Sec. IVlwas somewhat heuristic. 
We used ansatz (|63| without addressing the question whether it is self consistent, i.e., whether for VEVs B,C, D ^ 
the equation of motions will or will not lead to nonzero VEVs of other gluon fields. These questions will be addressed 
in the present section. The general analysis we will use is based on the symmetry consideration and the GL effective 
theory approach. 
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* T (iD M $*) + (h.c.)~4 8) 












TABLE IV: Composite vectors in the confinement picture. 



A. Symmetry breaking structure 



In this subsection, we consider possible symmetry breaking samples in the gauged NJL model and the structure 
of the corresponding homogeneous gluon condensates. Our strategy is the following. For each symmetry breaking 
sample, we will pick up the maximal set of gluon condensates consistent with it. The dynamics of course could allow 
a subset of this set to be a solution. This much harder issue is intimately connected with the structure of the GL 
effective theory corresponding to the symmetry sample and will be discussed in Subsec. IVI Bl 

In the present analysis, it will be convenient to use the (partly- unitary) gauge (JTTJ) for the diquark field Then 
the symmetry in the 2SC phase is that presented in Eq. f|43[) . In the general case, the homogeneous gluon condensates 

consist of 32 VEVs, (A^), (a = 1, 2, • • • , 8, n = 0, x, y, z). The symmetry (jj3} contains 9 parameters. However, since 
all gluon fields are singlet with respect to the chiral group [Z7(1) T 3 x U{l) T a] x , only 7 parameters connected with the 

group SU(2) C x f/(l) em x 50(3) ro t matter. By using the corresponding transformations, the 32 VEVs can be reduced 
to 25 ones. 

Let us show that these 25 VEVs can be chosen in the form: 

(4 8 >), (113a) 




and 




(Af), (£ — 1, 2, 3, j — x, y, z), (113c) 




(113b) 



(4 a) ) (a =1,2,..., 8), (113d) 



while (A { x s) ) = (A [ y s) ) = (4 5) ) = (4 7) ) = (4 4) ) = (Ai 5) ) = (£ ] ) = 0. 

The proof is going as follows. Because of the rotational symmetry, without loss of generality, we can choose one 
direction for the VEV of the spatial component of the 8th gluon, say, (A^) ^ 0. We may apply the SU(2) C symmetry 
to the same spatial component of 4>j and thereby obtain (4%) ~ (0 , (A^)). Note that we can still rotate <f> x and <p v 
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by a SO(2) spatial rotation around z-axis and vary the upper component of <fij by a U(l) transformation with the 
generator which is an appropriate linear combination of cr 3 in the color SU(2) C and the charge Q cm - By using these 
transformations, we can choose (<f> y ) to be real. As a result, we can reduce general homogeneous gluon condensates 
to the 25 VEVs indicated above. 

This set of 25 VEVs breaks symmetry P3"|) down to the chiral [U(1) T 3 x U(1) T ^] X [henceforth this irrelevant for our 
discussion chiral group will be omitted]. When we take subsets of the set in Eq. (|113| . typical symmetry breaking 
patterns are: 7 

SU(2) C x C/(l) em x SO(3) rot 

-» SU(2) C x U(l) em x SO(2) rot , (phase A), 

-» f/(l) cm x SO(2) rot , (phase B), 

-» SO(2) Mt , (phase C), 

-► S*0(2)dia g , (phase D), 

-»■ U{l) em , (phase E), 

— > nothing, (phase F). 

The maximal subsets of the gluon condensates consistent with these symmetry breaking patterns are 
(phase A) (4 8) )^0, (A^) + 0, 

(phase B) (4 3) ) ^ 0, (A { ^) f 0, (4 ?) ) ^ 0, (^ 8) ) f 0, (A< 3) ) f 0, (4 6) ) ^ 0, {Af) f 0, 
(phase C) (A { Q a) ) ± 0, (a = 1, 2, • • • , 8), (A«) ^ , (A®) ^ 0, <4 6 >) ^ 0, (A^) ± 0, 
(phase D) (4 2) >^0,(4 8) >^0, (^f) = (4 6) > ^ 0, 

(phase E) (4 3) ) ^ 0, (A { 6) ) f 0, (^ 7) ) + 0, (A< 8) ) ^ 0, (Af) f 0, (^f ) ^ 0, (A^) f 0, (4 8) ) + 0, 
(phase F) All 25 VEVs in Eqs. (lla)-(lld) ^ 

As was pointed in [l8|,[33[, the phase A corresponds to the single plane- wave LOFF phase [26]. In this case, both 
color electric and magnetic field strengths equal zero. This simplest case of the LOFF phase has been analyzed by 
several authors [26l. I28L f33(] . With the neutrality conditions taken into account, it was shown that the single plane- wave 
LOFF phase cannot resolve the chromomagnetic instability in the 2SC and g2SC regions [331 ] . 

While in the phase B the symmetry breakdown sample is the same as in the gluonic cylindrical phase II considered 
in Subsec. IV C[ in the phase C, it coincides with that in the gluonic cylindrical phase I discussed in Sec. [V] The 
VEVs of gluon fields in the latter constitute a subset of the VEVs presented in the maximal set in the phase C. Based 
on the GL approach, it will be shown in Subsec. IVI Bl that this subset is self consistent. On the other hand, as was 
already pointed out in Subsec. IV C[ the question concerning the existence of the gluonic cylindrical phase II is open. 
Recall that there are nonzero chromoelectric field strengths in these phases. 

In the phase D, the VEV set does not change under the SO(2)^ a „ transformations with the generator which is an 
appropriate linear combination of a 2 in the color SU (2) c and the generator of spatial y-z rotations. We may call this 



7 It is mathematically possible to consider also the phase with (j4q ) 7^ 0, (j4q ) 7^ and (A z ) 7^ 0, (A z ) 7^ 0, in which the unbroken 
subgroup is U(l) c x U(l) C m x 50(2) ro t with the generator a 3 for the color U(l) c - However, the Debye and Meissner screening masses 
for the SU(2) C gluons Aa ' do not take imaginary values both for 8fi < A and Sfi > A. Therefore there is no reason to consider such 
a phase. 
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(2) (8) 

phase the gluonic color-spin locked (CSL) phase. In this phase, while nonzero (A ) and (A ) are allowed, nonzero 

(Aq) and (Aq 3 ' : '' T ' > ) are not, because they break the S f O(2)di a g symmetry. It is noticeable that there exist both 
chromoelectric and chromomagnetic field strengths in the gluonic CSL phase. 

In the phase E, the U(l) em is unbroken, while the rotational 50(3) ro t symmetry is completely broken down. At last, 
the phase F contains the maximum number, 25, of gluon condensates. This classification provides a useful framework 
for studying dynamics of the gluon condensation. 

In conclusion, we comment on the relation between the homogeneous gluon condensates and inhomogeneous diquark 
condensates. It is clear that when there are no field strengths, constant gluon condensates can be removed by using 
an appropriate gauge transformation. However, the point is that such gauge transformations can break constraint 
(fTTj) for the diquark field. It is exactly what happens in the case of the single plane- wave LOFF state (the phase A): 
One can remove (A z ) at the cost of introducing an exponential factor with a linearly depending on z phase in the 
diquark field A. For other examples, see Ref. [341 ] . 

B. Ginzburg-Landau approach 

In this section, we study the GL approach for dynamics with vector gluon condensates in the gauged NJL model. 
As was shown in Sec. IIII1 in the hard dense loop approximation, the potential V in the model is given in Eq. (|32p . 
However, it is very hard to perform explicitly the calculations of the term with the fermion trace in V in the case of 
nonzero gluon condensates (|113|) . 

The GL approach can help to overcome this difficulty in studies of nearcritical dynamics in systems with a second 
order phase transition. In this approach, we need to pick up those local operators O n with the mass dimension four 
and less which are invariant under the symmetry of a system. In our case, it is 

SU(2) C x U(l) cm x [U(1) t b x C/(1) T |] X x 50(3) rot (114) 

(see Eq. ((12])). The GL effective action and Lagrangian, 5 e ff and £ c g, are written in terms of the sum of these 
operators: 

S oS [A°;fi,fjL e ,A] = /dx 4 £ cff [^; M ,/i e ,A], C cS = ^K n O n . (115) 

"* n 

The part of the Lagrangian without derivative terms yields the effective potential V e s. By analyzing the structure 
of the £ e ff, we can determine the relevant terms before providing the calculations of the full potential V (|32l) . It 
essentially reduces the amount of the work. 

Let us consider the GL effective Lagrangian in the gauged NJL model. Because the SU(2) C in Eq. pi4[) is a gauge 
symmetry, the building blocks of C c e are 

4> , fa, Do, Vj, fa, f jh , Af\ (116) 

where the indices j and k represent spatial components. The coefficients K n of the operators are functions of A, /j,, 
fx e , and fig. But when we expand the effective action with respect to fig, one should add the field A : 

4> , fa, Do, Vj, fa, f jk , 4 8) , Af\ (117) 

and now the the coefficients of the operators are functions of A, fi, and /i e . For completeness, we will describe the 
larger set of the operators expressed through building blocks (|117[) . 

Up to the mass-dimension four, we find the operators without the singlet field A^ , 

4>Ua, (118a) 
(f>l(Wo^o), <^(*W), (H8b) 

(40o) 2 , (4<M(^<n (4^M<h), w>y') 2 , (</>^ fe x4<n o^x^), (use) 
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|2? </>o| 2 , |P^ | 2 , l^l 2 , l^^fcl 2 , ^((A^ + ^Z^'.Z)*})^, 4({^> }^), (H8d) 

and 

tr(/ 0j /° J ), te(fikf ik ), i<t>)P k <t>k, <t>\f j <t>i- (H8e) 

(8) 

Including , we also get the following operators, 

^0o4 8) - $^'4 8) . ^^o(4 8) ) 2 - ^}^(4 8) ) a . ^o^a 4 8) , $4? hA!®, 4^d 3 Af\ (119a) 

0S(ii?o0o)4 8) . 0}(W)4 8) , 4(^^)4 8) , (nob) 

and 

0j^4 8) , 4^4 8) 4 8) , i&M^) 3 , ^(4 8) ) a , jt /4 8) ^i 8) ' <f>WdoA?\ <t>]<t> k VAf, (120a) 

4<hd*A?\ ^d k Af\ Afm(iV Q( jj), Af4{iV^ Q ), Af^](iV k ^), Af^{iV^ k ). (120b) 
The operators including only the singlet field are 

4 8 \ (4 8) ) 2 , (4 8) ) 2 = 4 8) (4 8) ) 2 , ete -- ( 121 ) 

The effective Lagrangian is given by the sum of the above operators. It is understood that the hermitian conjugate 
operators are also included in £ c s, and some coefficients can be complex. The following remarks are in order, a) As 
has to be, the sets (|118p . (|119p . and (|120|) do not contain operators corresponding to the Debye and Meissner screening 
mass terms for the SU(2) C gauge gluons. b)We do not consider the parity violating term such as ^ k </>lfoj^>k , etc.. 
c) It is easy to check in the hard dense loop approximation (with the one loop fermion contributions) that due to 
the structure of the fermion-antifermion-gluon vertex, a n-point vertex of gluons has an even number of the e b -type 
vertices. Therefore we do not consider terms as e Q/3 ((/>j) Q (i2?o < / , ' : ')/3 ,etc, where a,f3 denote the SU(2) C indices. 

We are ready to show that the GL effective potential has the form presented in Eq. ([64]) . Indeed, now we know all 
relevant and marginal operators which can be constructed for the fields B,C and D defined in Eq. They are: 

<j>)<fp corresponding to B 2 , $\{iV a &) to DB 2 , \V 3 (j> k \ 2 to B 2 C 2 , |£> ^| 2 to B 2 D 2 , tr{f 0j f 0j ) to C 2 D 2 , and (^') 2 , 
|(/)J^ fe | 2 and (4>j(j> k ) 2 corresponding to B 4 . Since there are no other relevant and marginal operators which can be 
constructed from these three fields, we are led to the form in Eq. f|(34[) . 

To obtain all of the coefficients K n of the operators is a very hard task. Only some of them are known. The coeffi- 
cients of <Pq4>o , (4^) 2 j an d ^j'P j (4 ) 2 corres P on d to the Debye and Meissner screening masses, respectively [l3| . 
The coefficients of the operators ^(iDq^P), \T>j(j)k\ 2 , and tr(/(y/ 0j ) correspond to the parameters Tub, Xbc, an d 
^cd, respectively, in the GL potential in the gluonic phase used in Sec. [V] (see Eq. (fTo]0 . These parameters are 
calculated in Appendix [B] Without 5/j,, the three and four point gluon vertices have been obtained in Ref. [48|. The 
important consequences of those calculations are that the coefficients of the three point vertex operators <pQ(iT>Q(j) ), 
(/>](i2?o0 : ')j arid 4> (iT>j(j>!) are vanishing at the order of /i 2 , (i.e., they ~ O(fi)), and that four point vertex operators 

such as \T>j<j)k\ 2 , tr(/oj/° : '), and A^ (j)\{iD k 4P) are of order fi 2 . For their concrete expressions, see Ref. (48|. The 
results of our calculations in Appendix IB] are of course consistent with them. 

The GL approach can be applied only to those symmetry breaking samples, discussed in Subsec. I VI Al which are 
connected with a second order phase transition. Therefore, it cannot be applied both to the sample A (the single 
plane- wave LOFF state) and the sample B (the gluonic cylindrical phase II) , which are related to a first order phase 
transition (see Refs. (2a. [28l [33l] and the discussion in Subsec. IVI A[) . 8 On the other hand, as was shown in Sec. [Vj 



Note that as follows from Refs. |l3t Il7| . the instabilities at A = 8fi in the g2SC phase (not considered in this paper), also cannot be 
connected with a conventional second order phase transition. 
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the GL approach consistently describes the gluonic cylindrical phase I, assigned to the symmetry breaking sample C 
(a detailed analysis leading to this conclusion is considered in this subsection below) . The discussion of the dynamics 
in the symmetry breaking sample D (the gluonic CSL phase) is beyond the scope of this paper. 

To demonstrate the power of the GL approach, we will apply it to prove a self-consistency of the ansatz for the 
gluonic phase we used in Sec. [V] In Subsec. IV CI a tadpole mechanism for producing VEVs of gluon fields from triple 
vertices was discussed. Let us describe it in more detail. One can always divide a field into the VEV part and the 
fluctuation one, 

= (Af)+af, (£-1,2,3), (122a) 
0m = (0/»> + C/., (122b) 
Af = {Af)+af. (122c) 

Substituting expressions (|122|) for all fields into the GL potential expressed through operator sets (I118p - (112ip . we 
obtain the VEV part, the tadpole part, which is linear with respect to fluctuations, plus higher order terms. The 
gap equation is the stationary condition for the effective potential and it is equivalent to vanishing the tadpole part. 
We emphasize that this condition is stronger than solving the gap equation for a given ansatz by taking into the 
consideration only the fluctuations of the fields in the ansatz and ignoring the fluctuations of other ones. 

The condition of vanishing the tadpole contribution has to be correct for any consistent ansatz used for VEVs: its 
violation indicates that the ansatz is inconsistent. While for the maximal ansatz, corresponding to a chosen symmetry 
breaking pattern, the violation would imply that the pattern cannot be realized in the model at all, for a non-maximal 
one, this could just mean that the latter is not closed and one should add new VEVs for restoring its consistency. 

As a concrete example, let us consider the ansatz with (A^) ^ 0, (A^) ^ 0, and (Ay 8 * 1 ) ^ 0, which was discussed 
in Ref. [l^], and show that it is not self-consistent. The point is that the operators A^ cj^AiD^cj^) and ^{iDq^) 
yield the tadpole contributions a y 3 \A y 8 ^)(A^) 2 and (A^ ) 2 , respectively. Because there are no other tadpole 
terms for a y 3 ^ and the condition of vanishing the tadpole contribution leads to (A^) = 0, which is inconsistent 
with the original ansatz. If we started from the ansatz with (A^) ^ 0, (A^) ^ 0, (A y 3 ^) ^ 0, (A^) ^ 0, and 

(Ay ) 7^ 0, new tadpole contributions for a y and a would appear. They will either allow to cancel the previous 
tadpole contributions without generating new ones or this modified ansatz is also inconsistent and should be in its 
turn extended. This process will eventually lead either to a consistent ansatz or to the conclusion that this symmetry 
breaking pattern is not realized in the model. 9 

Without calculating the coefficients of the operators and solving the gap equations, this check helps to pick up a 
self-consistent ansatz. By estimating the size of the coefficients, one can also specify suppressed VEVs, as we will see 
below. The latter can be useful for simplifying the ansatz we work with. 

Let us now turn to the gluonic phase in Sec. [V] As was shown in Subsec. VI A, the phase C, having the same 
symmetry structure as the gluonic cylindrical phase I, can contain other VEVs than B = g(A < f ' > ), C = g(A^P), 

D = g(A^), and fi$ — \fi/2g(A^) used in Sec.V. Is the ansatz including only B, C, D, and [i$ self-consistent? 
Based on the GL approach, we will show that the answer to this question is affirmative. Actually, we will show that, 
strictly speaking, the self-consistent minimal ansatz for the gluonic cylindrical phase I is 

(A^),{A^),{A^),{A^), and (A { 4) ), (A^) . (123) 

However, it will be shown that the additional VEVs (A^) and (A^) are suppressed in the vicinity of the critical 
point. 

For B,C',D 7^ 0, the tadpole contributions for a[ ) 4:5 ' ) come from the operators (^(iVjcf^) and (f>of°^(f>j in the GL 
effective Lagrangian with, say, coefficients K\ and K2, respectively. The condition of vanishing the tadpole part yields 
D ~ —Ki/K2- On the other hand, the solution of the gap equations is given in Eq. ([71]) . Since M% and Tub are 
the coefficients of different operators, 4>^<j)' and (/^(Wo (/)•'), respectively, it is hard to expect a "magic" relation (fine 



This example corresponds to one of the non-maximal ansiitze for the symmetry breaking sample E considered in the previous subsection. 
For the choice of the VEVs in the phase E used in the present paper, {A^) 0, (A^) ^ and (Ay') ^ in the ansatz should be 
replaced by {Af ] ) ^ 0, (A ( y 6) ) ^ and (a{ 8) ) ^ 0. 
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tuning) among K\, K2, Mj, and Tub coefBcients, which would lead to the same solution for D from the gap equation. 
Therefore, one should introduce VEVs for A^' 5 ^ in this case. 

Noting that A 1 ^'' 5 ^ have nonzero Debye mass )tc,d 4 ~ O(n), one can then estimate the solution of the gap equations 

for (Aq 4,5 -*) as (Aq 4 -*), (A^) ~ K\BC /m D 4 + K2BCD /m D 4 . It is obvious that near the critical point, where 
B,C,D — > 0, the contribution of the first operator dominates. Without S/j,, the coefficients K\ and were calculated 
in hard dense loop approximation in Ref. [481 ]. It was shown there that while K 2 is of order /z 2 /A 2 , K\ is vanishing 
in this leading approximation. It implies that K\ ~ O(fi) or less, and thereby we conclude that the contribution of 
this operator to the effective potential is the term B 2 C 2 with the coefficient ~ 0(1) at most. On the other hand, the 
coefficient of the B 2 C 2 term discussed in Sec.V is Xbc in Eq. (|85|) : Xbc ~ Oip 2 1 A 2 ). Therefore, the contributions 
of (A^) and (A^) can be ignored in the vicinity of the critical point. In addition to this, one can show that with 
the ansatz (I123[) for VEVs, there are no potentially dangerous operators in Eqs. (|118[) - (|12ip which would lead to 
tadpole contributions for A\ and other additional fields from the maximal ansatz absent in Eq. (|123[) . This allows 
us to conclude that the ansatz used in Sec. [V] is self-consistent indeed. The above example clearly shows the power 
of the GL approach and its relevance for the dynamics with vector condensates of gluons. 10 



VII. CONCLUSION 



The gluonic phases yield an example of dynamics in gauge models with matter in which the Higgs mechanism 
is provided by condensates of gauge (or gauge plus scalar) fields. Because VEVs of spatial components of vector 
fields break the rotational symmetry, it is natural to have a spontaneous breakdown both of external and internal 
symmetries in this case. Dynamics in such systems are quite sophisticated. The existence of exotic hadrons in the 
gluonic phases is especially intriguing. 

What could be directions for future studies in these phases? It is evident that it would be interesting to consider 
the spectrum of light collective excitations there. The results in Refs. [H,|4l[, obtained in the gauged a- model with a 
chemical potential for hypercharge (briefly discussed in Sec. Hi]) , suggest that the spectrum should be rich, containing, 
in particular, gapless NG modes, rotonlike and vortexlike excitations. 

Another interesting direction would be to clarify whether the symmetry breaking patterns considered in Subsec. 
IVI A[ can be realized as stable, metastable or unstable ground states in dense QCD. Recent results [49| . showing 
that the landscape of such ground states in the gauged cr-model with a chemical potential for hypercharge is rich, are 
encouraging. 

It is clear that it would be worth to figure out whether phases with vector condensates of gluons could exist in 
dense matter with three quark flavors. Recently, this possibility has been mentioned in Refs. [lTl [l9l. l50l. l5~fl] . 

It has been recently revealed that a strong enough magnetic field can influence the phase structure in dense quark 
matter (52|. 11 It would be interes ting to study this phenomenon in gluonic phases, especially because one can expect 
the existence of vortices there [U l55| 

Last but not least, it would be worth searching for a realization of a gluoniclike phase in condensed matter. Recently, 
there has been a considerable interest in systems with coexisting order parameters (such as high T c superconductors) 
in condensed matter (56j. Generating vector condensates is a natural way of creating such systems (for example, in 
the gluonic cylindrical phase I, electric superconductivity coexists with spontaneous rotational symmetry breaking). 
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10 Another example of this analysis is for the gluonic CSL phase D. One can show that with the GL operators in Eqs. (118)-(121), there 
are two self consistent ansatze in that case: the maximal one and the maximal one minus the VEV of the field A . 

11 This phenomenon has some similarities with the phenomenon of the magnetic catalysis in vacuum field theories |5.'j|. in particular, in 
the vacuum QCD [53 ]. 
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APPENDIX A: STRUCTURE OF NEARCRITICAL SOLUTION IN GLUONIC PHASE 



In this Appendix, we analyze the structure of the nearcritical solution with B =/= 0, C =/= 0, and D =/= discussed in 
Subsec. [V3 Eqs. §^-§7} yield 
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On the other hand, Eq. (|68|) yields: 
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Because the coefficients M B , Tdb, etc. in Eq. (|A1[) are the functions of A, /z e and /x 8 , Eqs. (|A1|) and (|68|) constitute 
a coupled system of six equations. However, as will be shown below, near the critical point, the equations (|A1|1 and 
(|68|) decouple. 

Near the critical point, M B ss 0, one can expand A, /i e and /i 8 around B = C = D = 0: 
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(A3) 
(A4) 
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where the bar-quantities correspond to the 2SC solution with B = C = D = 0. By definition of the bar-quantities, 
the following stationary conditions are satisfied: 
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A = A,/i E =/i e ,^8=^8 



They yield 
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Note that the bar-quantities A, /2 e and p,g are uniquely determined by the theoretical parameters G A ,A,[i in the 
gauged NJL model. Expanding the left hand side in Eq. (|A2[) around the bar-quantities, we obtain the gap equation 
for £a, £e and £ 8 expressed in a matrix form as 
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where all derivatives are calculated at A = A, /i e = fx e , and fig = /ig. 

Combining Eq. (|A"Tj) with Eq. (|A"TU|) . we find the following approximate solution: 



sol 



-Ml 



-A 



CD 



3|Tdb| y 



-Ml 



3A 



sol 



-Ml 



£Sol £Sol £Sol 

i ?8 



£SOl £S 

S"A )Se 



0{{M B )% 



(All) 



where higher order terms in Ml were neglected. Here all the coefficients M B) Tdb, Xbc, and Xcd are expressed 
through the 2SC values A, /2 e and fig. This solution exists when 



Ml < 0, A BC > 0, Xcd < 0. 



(A12) 



Note that in Eq. (|Af 1[) the convention B so i > and C so i > is chosen. 

Substituting solution (|A11|) into GL potential (|64|) , we find that its gluonic part V e g — V A consists of leading terms 
~ C((Ml) 3 ) and subleading terms ~ 0((M B ) ). In particular, the deviation of Va from that in the 2SC solution is 
of subleading order 
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Indeed, the second order of the Taylor expansion for Va yields 
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where Eq. (|A10j) was used (here all derivatives are calculated at A = A, fi e — /i e , and /is = fis)- 
Omitting the subleading terms, we arrive at the reduced effective potential, 



Veft = V A (A,fl e ,fl 8 ) + -M 2 B 2 +T DB DB 2 



XbcB 2 C 2 



XcdC 2 D 2 . 



(A15) 



This potential is composed of two parts: the "constant" 2SC part Va, with frozen fermion parameters, and the 
dynamical gluonic part. 



APPENDIX B: COEFFICIENTS IN GL EFFECTIVE POTENTIAL 



In this Appendix, we calculate the coefficients Xbc, Xcd, and Tdb of the marginal and relevant operators in 
the reduced GL effective potential ([78]) . The coefficients Xbc and Xcd are connected with the operators \Vj4>k\ 2 
and ti^/iy/ - 7 '), respectively. Therefore, they can be obtained from the kinetic terms of the 4- 7th and l-3rd gluons. 
Although the coefficient Tub of the triple vertex is connected with the operator (j)UiT>o^), it comes from a non- 
hard-dense-loop part, as we will see below. Because of that, Tdb will be directly calculated from the corresponding 
three-point vertex. In passing, the coefficient Xb in GL potential is connected with the operator ((pj^) 2 and 
cannot be reduced to a calculation of a two-point function. On the other hand, because the coefficient Xb d in (|64|) is 
connected with the operator \T>Q(j)j\ 2 , it can be obtained by taking time-derivatives of the vacuum polarization tensor 
for the 6th gluon. 
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One of the main difficulties in calculating these coefficients is related to the presence of the gauge fields C = gA^P 
and D = gA^ connected with noncommutating color matrices T 1 and T 3 . In order to overcome it, the following 
trick will be used. We decompose the inverse Nambu-Gor'kov propagator (f2T)|) with gluon fields as 



where 



S- 1 (X) = Sc 1 (X)+Mbd, (Bf) 



and 



/ DT 3 7° - BTV 
I -L>T 3 7 ° + 5TV 



with X^ = (x , a; , x 2 , x 3 ) = (i, x, y, z) denoting space-time coordinates. The point is that since the diquark gap A is 
S^(2) c -invariant, we can remove the field C = gA% from Sc(X) by using a SU(2) C gauge transformation. [Note that 
since the field strength is nonzero in the gluonic phase, it is impossible to remove all gauge fields from the propagator 
S g (X)]. In fact, one can easily find this SU(2) C gauge transformation: 

/ -iCT^x 3 Q \ 

U W S [ e^ 3 ) (B4) 
transforms the inverse S^ 1 (X) in Eq. (|B2[) into 5 f_1 (X) without C, i.e., 

Sc 1 (X) = U*(X)S- 1 (X)U(X). (B5) 
Let us now expand the potential V (]29[) in the power series with respect to Sc{P) in momentum space: 

00 (—'\) n f d 4 P 

V = V A + V g+ ^2^J-j —-T^SoM^r, (B6) 

where the identity 

Tr In Sc 1 = Tr m( US^rf ) = Tr In 5 _1 (B7) 

was used. Note that, by construction, the V g term is C independent. 

In order to utilize expression (|B6[) . we need to calculate the propagator Sc{P) in momentum space. From Eq. (|B5|) 
we get 

S C (X) = U*(X)S(X)U(X). (B8) 



Using then the relations 
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a ho, (bio) 



we find that the Fourier transform of e lC a x is given by the sum of two 5-functions: 

FTe iC'^x* = 1 ( §(p 3 + a) + S{p 3 _c')) + V ( 5{p 3 + C) - 5{p 3 - C) ) . (Bll) 
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This implies that the Fourier transform of U(X) (|B4[) is expressed through these 5-functions. And since the Fourier 
transform of S(X) is known, the Fourier transform of Sc(X) can be easily found from Eqs. (|B8[) and (jBlip : 
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While G^ is given in Eq. (|21|) . the matrix elements G^„ and in S rg {P) are connected with the gapped part S A (P), 
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where the quantities G A , G^, Si2,2i are given in Eqs. (|20|) . (23), and (|24|) . In fact, G^r g and are: 
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where Pc is the four-vector (0,0, 0,G/2). Note that S^ct 1 = — cr 1 E! A and, using Eq. (Tj"5]) , one can check that 

S-HP)S rg (P) = l. 

After taking the trace over the color indices, we obtain the square term for Xbc and Xcd , 
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In fact, the B 2 - and D 2 -terms in Eq. (|B19|) correspond to the vacuum polarization tensors II44 55 and II™, respectively 
(see Eqs. (76) and (40) in the second paper in Ref. [HI-) Therefore, the calculations of Xbc an d Xcd are reduced to 
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those of U J 4 J 4 55 and 11$. We find 
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In order to calculate the coefficient T^b, the cubic term in expansion (|B6[) is required: 
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The functions IC eie2e3 and Ti. eie2£3 are expressed through the following components of the propagators (compare with 
Eqs. (20)-(24)): 
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where r = ±1 are eigenvalues of T3, and 

H±(P) = -iA[3± r= _ 75 A++H± T= _ 75 A;], 
H±(P) = -iA[3± r=+7s A+ + H± T=+75 A;] 



(B26) 
(B27) 



(B28) 
(B29) 
(B30) 
(B31) 



(B32) 
(B33) 



with 



±\2 ' 



The explicit expressions of K. eie2e3 and W eie2e3 are: 
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where we used the relations 
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After integrating over p°, we find 
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The integrand I2 contains only the particle-antiparticle contribution and hence should be negligible. We expand ii 2 
in the series with respect to /j 8 and obtain 
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As was expected, the contribution of I2 is suppressed indeed 
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